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Abstract 

It is fairly well known that rotation in three dimensions can be expressed as a quadratic 
in a skew symmetric matrix via the Euler-Rodrigues formula. A generalized Euler-Rodrigues 
polynomial of degree 2n in a skew symmetric generating matrix is derived for the rotation 
matrix of tensors of order n. The Euler-Rodrigues formula for rigid body rotation is recovered 
by n = 1. A Cayley form of the n*^ order rotation tensor is also derived. The representations 
simplify if there exists some underlying symmetry, as is the case for elasticity tensors such 
as strain and the fourth order tensor of elastic moduli. A new formula is presented for the 
transformation of elastic moduli under rotation: as a 21-vector with a rotation matrix given by 
a polynomial of degree 8. Explicit spectral representations are constructed from three vectors: 
the axis of rotation and two orthogonal bivectors. The tensor rotation formulae are related to 
Cartan decomposition of elastic moduli and projection onto hexagonal symmetry. 

1 Introduction 

Rigid body rotation about an axis p, |p| = 1, is described by the well known Euler-Rodrigues 
formula for the rotation matrix as a quadratic in the skew symmetric matrix P, 



^ ' Pij = -eijkPk, (1) 

« ! 

Q Where eijk is the third order isotropic alternating tensoi0. Thus, 

Q = exp(^P) =I + sin^P + (l-cos^)P^ (2) 

X : 



where 6 is the angle of rotation. Euler first derived this formula although Rodrigues (1211 ) obtained 
the formula for the composition of successive finite rotations (fill : 0). The concise form of the 
Euler-Rodrigues formula is basically a consequence of the property 

p3 = -P. (3) 



^Thc summation convention on repeated subscripts is assumed in ([T])i and most elsewhere. 
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Hence, every term in the power series expansion of exp(^P) is either P^'^ = (— 1)'^'"'"^P^ or p^'^+i = 
(— 1)'^P, and regrouping yields ([2]). Rotation about an axis is defined by the proper orthogonal rigid 
body rotation matrix Q(6', p) G 5*0(3), the special orthogonal group of matrices satisfying 

QQ* = Q*Q = I, (4) 

and the special property det Q = 1. The generating matrix P(p) G so(3) where so{n) denotes 
the space of skew symmetric matrices or tensors, usually associated with the Lie algebra of the 
infinitesimal (or Lie) transformation group defined by rotations in SO{n). Alternatively, 5*0 (3) is 
isomorphic to so(3) via the Cayley form (jj) 

Q = (I + S) (I - S)-' = (I - S)-i (I + S), (5) 

where S G so(3) follows from ([2]) as 

S = tanh(f P) = tanf P . (6) 

The first identity is a simple restatement of the definition in ([5]) with Q = exp(6'P), while the 
second is a more subtle relation that depends upon ([3]) and the skew symmetry of P, and as such 
is analogous to the Euler-Rodrigues formula. 

We are concerned with generalizing these fundamental formulas for rigid body rotation to the 
rotation matrices associated with Q but which act on second and higher order tensors in the 
same way that Q transforms a vector as v — ^> v' = Q*v. If one takes the view that the rotation 
transforms a set of orthonormal basis vectors 61,62,63, then the coordinates of a fixed vector 
V = viei + 1>262 + ^363 relative to this basis transform according to {vi,V2, Vs) {v[,V2, v'^), where 
= QijVj. A tensor T of order (or rank) n has elements, or components relative to a vector basis, 
Tij,,,ki- Under rotation of the basis vectors by Q, the elements are transformed by the relations 

'^ij...kl ~ Qij...klpq...rsTpq,,,rsj (7) 

where Q = Q(^, p) is a tensor of order 2n formed from n combinations of the underlying rotation, 

^ij ■ ■ ■ kl pQ . . . rs ~ QipQjg ■ ■ ■ QkrQis • (8) 



Q is also known as the n-th Kronecker power of Q (1l8l). Lu and Papadopoulos (jl4l ) derived a 
generalized n-th order Euler-Rodrigues formula for Q as a polynomial of degree 2n in the tensor 
P, also of order 2n, defined by the sum of n terms 

Pij...klpq...rs ^ip^jg " " " ^kr^ls ~l~ ^ipPjq ■ ■ ■ ^kr^ls ~l~ • • • ~l~ ^ip^jq ■ ■ ■ Pkr^ls ~l~ ^ip^jq ■ ■ ■ ^krPls • (9) 

Thus, 



exp(^P) = P2n(^, P), (10) 



where 



n 



P2„,(^, x)= V e^^^ n (^^), withAfc = A:^, i = ^. (11) 
, M — -^j 

k=—n j=—n ■' 

Lu and Papadopoulos' derivation of ffTTl) is summarized in Section [31 The case n = 2 has also been 



derived in quite different ways by Podio-Guidugli and Virga ( l20l ) and by Mehrabadi et al. (jl6l ). 



^Q* denotes transpose and I is the identity. 
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The purpose of this paper is twofold. First, to obtain alternative forms of the generalized 
Euler-Rodrigues formula and of the associated Cayley form. Using general properties of skew 
symmetric matrices we will show that the n-th order Euler-Rodrigues formula can be expressed in 
ways which clearly generalize the classical n = 1 case. The difference in approach from that of Lu 



and Papadopoulos (jl4l ) can be appreciated by noting that their formula (fTTj) follows from the fact 
that the eigenvalues of P are A^, k = — n, . . . ,n, of eq. (fTTI) . The present formulation is based on 
the more fundamental property that P has canonical form 

P= Pi + 2P2 + ... + nP„, (12) 

where Pj satisfy identities similar to ([3]), and are mutually orthogonal in the sense of tensor product. 
This allows us to express exp(^P) and tanh(|P) in forms that more clearly generalize ([2]) and ([6]). 
Furthermore, the base tensors Fj can be expressed in terms of powers of P, which imply explicit 
equations for Q and the related S. We will derive these equations in Section [31 

The second objective is to derive related results for rotation of tensors in elasticity. The under- 
lying symmetry of the physical tensors, strain, elastic moduli, etc., allow further simplification from 
the 2n— tensors Q G S'0(3"') and P G so(3"'), where the nominal dimensionality 3" assumes no 



symmetry. Thus, Mehrabadi et al. (jl6l ) expressed the rotation tensor associated with strain, n = 2, 
in term of second order tensors in 6-dimensions, significantly reducing the number of components 
required. By analogy, the rotation matrix for fourth order elastic stiffness tensors is derived here 
as a second order tensor in 21-dimensions Q(6',p) = exp [OP) G 5*0(21) where P G so(21). This 
provides an alternative to existing methods for calculating elastic moduli under a change of basis, 
e.g. using Bond transformation matrices (0) or other methods based on representations of the mod- 
uli as elements of 6 x 6 symmetric matrices ffl6l ). Viewing the elastic moduli as a 21- vector is the 
simplest approach for some purposes, such as projection onto particular symmetries (flil: IgI). In fact, 
we will see that the generalized Euler-Rodrigues formula leads to a natural method for projecting 
onto hexagonal symmetry defined by the axis p, as an explicit expression for the projector appears 
quite naturally. 

The third and final result is a general formulation of the so-called Cartan decomposition (0) of 
tensor rotation. The action of Q of ([2]) on a vector leaves the component parallel to the axis p 
unchanged, and the component perpendicular to it rotates through angle 6. The latter decreases 
the part perpendicular to p to cos^ of its original value, and introduces sin 6* times the same 
magnitude but in in a direction orthogonal to both the vector and the axis. This simple geometrical 
interpretation generalizes for tensors, in that we can identify "components" that remain unchanged, 
and others that rotate according to cos j9, sin j9, j = 1,2, ... ,n. The components form subspaces 
that rotate independently of one another, just as the axial and perpendicular components of the 
vector do for n = 1. This is the essence of Cartan decomposition of tensors, see (jsj) for further 
details. Our main result is that the components and the subspaces can be easily identified and 
defined using the properties of P and Q that follow from (fT2l) . 

We begin with a summary of the main findings in Section [2], with general results applicable to 
tensors of arbitrary order described in Section [31 Further details of the proofs are given in Section 
[H Applications to elasticity are described in Section [5l where the 21 dimensional rotation of elastic 
moduli is derived. Finally, the relation between the generalized Euler-Rodrigues formula and the 
Cartan decomposition of elasticity tensors is discussed in Section [61 

A note on notation: tensors of order n are denoted P, Q, etc. Vectors and matrices in 3D are 
normally lower and capital boldface, e.g. p, Q, while quantities in 6-dimensions are denoted e.g. 
Q, and in 21-dimensions as Q. 
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2 Summary of results 

Our principal result is 

Theorem 1 The rotation tensor Q(^, p) = exp (^P) has the form 



Q= 1 + ^ (sinA:^Pfc + (l-cosA;e)P^) (13) 
fe=i 

where P^, k = 1,2,. ..,n are mutually orthogonal skew-symmetric tensors that partition P and have 
the same basic property as P in 

P= Pi + 2P2 + ... + nP„, (14a) 
PiP,- = P,Pi = 0, z^j, (14b) 
P,^ = - Pi . (14c) 

The Cayley form for Q is 

Q= (1+ S)(I- S)-^ = (I- S)"^(I+ S), (15) 

where 

n 

§ = ^ tan(A;f ) P^ . (16) 

k=l 

The tensors P^ are uniquely defined by polynomials of P of degree 2n — 1, 

" 2 I -2 

Pfc = Pn,fc(P), where PnA^) = j: YlijTli^) ^ 1 < k < n. (17) 

The polynomial representation Q = P2n(^) P) has several alternative forms, three of which are as 
follows: 

" 2('— " 

P2n(^, x) = 1 + V Vttt TTT sin A;^ x + (1 - cos A;^) x'l FT (x' + j') (18a) 

^ (n-/c)!(n + A;)! ^ 

n /r) • 9N2fc— 1 fc— 1 



where 



l + sinMl + ^(l + f + + 
+ (1 - cos^)x^(l + ^(1 + ^(1 + ^(1 + . . . ^)))) , (18c) 



X, = (1 - cos^)(x2 + j2)/(2j + 1). (19) 



We note that n = 1,2,3,4 corresponds to rotation of vectors, second, third and fourth order 
tensors, respectively. All these are covered by considering n = 4 in Theorem [T], for which the series 
fll8cp and fll8ap are, respectively. 



Ps{e,x) =l + sin^x( l + ^^^^^(x2 + l) 



1 — cos6'), r. , (1 — COS 



X 



+ 4)(l + ^^-^(x^ + 9)) 



2.5 ' 3.7 
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I i — cosy I, 2 A\/, (i — cosy I, 2 
1 + 3^ ^ (a:^ + 4) (1 + ^ (o:^ + 9)) 



= 1 + [ sin ^ X + (1 - cos (x^ + 4) (x^ + 9) (x^ + 16) 

- ^ [2 sin 2^ X + (1 - cos 2^)x21 (x^ + 1) (x^ + 9) (x^ + 16) 
6! 

+ ;^[3sin3ex + (1 -cos3^)x2](x2 + l){x'^ + 4){x'^ + 16) 

- ^ [4 sin 4^ X + (1 - cos 4:9)x'^] (x^ + 1) (x^ + 4) (x^ + 9) . (20) 

This expression includes all the Euler-Rodrigues type formulas for tensors of order n < 4, on account 
of the characteristic polynomial equation for P of degree 2n + 1, 

P(p2+ I)(p2 + 4I)...(p2 + n2l) =0. (21) 

Thus, the last line in (l20l) vanishes for n = 3, the last two for n = 2, and all but the first line for 
n = 1, while the terms that remain can be simplified using (12T]) . The formulas for P2n, n = 2 and 
3, are 

PJ9,x) =1 + sm9x(l + -(1 - cos^)(x^ + 1)) + (1 - cos^)x2(l + -(1 - cos^)(x2 + 1)), (22a) 

3 6 

Pe{9,x) =1 + sin^x(l + -(1 - cos9){x^ + 1) [l + — (1 - cos^)(x2 + 4)]) 

3 10 

+ (1 - cos^) x^(l + -{1 - cos9){x^ + 1) [1 + —(1 - cos^)(x2 + 4)1) . (22b) 
6 15 

These identities apply for n < 2 and n < 3, respectively, since they reduce to, e.g. the n = 1 
formula using (12T1) . 

The second main result is a matrix representation for simplified forms of P and Q for 4th order 
elasticity tensors. Although of fourth order, symmetries reduce the maximum number of elements 
from 3^ to 21, and hence the moduli can be described by a 21-vector and P and Q by 21 x 21 



matrices. This reduction in matrix size is similar to the 6x6 representation of Mehrabadi et al. fll6l ) 
for symmetric second order tensors. The case of third order tensors (n = 3) is not discussed here 
to the same degree of detail as for second and fourth order, but it could be considered in the same 
manner. The following Theorem provides the matrices for n = 1,2,4, corresponding to rotation of 
vectors, symmetric second order tensors, and fourth order elasticity tensors, respectively. 

Theorem 2 The skew- symmetric matrices and the associated rotation matrices defined by the ten- 
sors (P, Q) for n = 1,2,4, are given by (P,Q), (P, Q) and (P, Q), respectively. The m x m, 
m = 3, 6, 21, skew symmetric generator matrix for each case has the form 

P = R-R\ (23) 

where 

R = — X, for vectors, n = 1, (24a) 

R = f ^2 ^'^^1 ' ■^^^ symmetric tensors, n = 2, (24b) 
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tensors, n = 4, 



(24c) 



with = Oaxs and 



X 



P3 0' 
Pi 

yP2 



P2 0' 
P3 

<Pi 



■ Pi 

P2 
<P3 



N 



^pi 

P2 

ps, 



(25) 



/n each case, the Q-matrix is given by Theorem [21 and the tensors acted on by Q are m-vectors 
which transform like v' = Qv. 

The specific form of the 6- and 21-vectors are given in Section [5l 

The third and final result identifies subspaces that are closed under rotation: 



n. 



Theorem 3 The Cartan components of an n-th order tensor T are defined as 

To = Mo T, Tj = Mj T, = Fj T, j = 1, 2 
where the n + 1 symmetric projection tensors, M^, k = 0,1, . . . ,n, are 

Mo = I + Pi + P2 + . . . + P^, Mi = - Fl 
They satisfy, for < k < n and 1 < i < n, 

i 7^ k. 



1,2, 



n. 



life = iVllfc . 



k= MfeP, = PiMfe = 0, 
M,Pi = P,Mi = Fi, 

i-ir (2 - 4o) 



(n - k)\{n + k)\ 



(26) 



(27) 



(28a) 
(28b) 

(28c) 



The projection Tq along with the n— pairs { Tj, Mj}, j = 1, 2, . . . , n define n + 1 subspaces that 
are closed under rotation: 



>Tn = Tn 



iTj = COS Tj + sin j^Mj, 



cosj9Rj-smj9Tj, j = l,2,...,n. (29) 



The meaning will become more apparent by example, as we consider the various cases of tensors of 
order n = 1,2, and 4 in Section O 
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3 General theory for tensors 
3.1 P's and Q's 

The transpose of a 2n-th order tensor is defined by interchanging the first and last n indices. In 
particular, the skew symmetry of P and the definition of P in ([9]) implies that P* = — P. Based on 
the definitions of Q in ([8]) and the properties of the fundamental rotation Q G SO (3), Q of (fTOl) 
satisfies 

^ = PQ, (30) 

where the product of two tensors of order 2n is another tensor of order 2n defined by contracting 
over n indices, {AB)ij„Mpq...rs = Aj...kiab...cdBab...cdpq...rs- This gives meaning to the representation 

Q = e'^. (31) 

Based on the skew symmetry of P it follows that 



= I, where Iij...klpq...rs = ^ip^jq ■ ■ ■ Skr^ls ■ (32) 

The isomorphism between the space of n-th order tensors and a vector space of dimension 3^^ implies 
that Q e SO (3") and P e so(3"). 

The derivation of the main result in Theorem [1] is outlined next using a series of Lemmas. Details 
of the proof are provided below and in the Appendix. 

Lemma 1 For any N > 1, a given non-zero N x N skew symmetric matrix B has 2m < N 
distinct non-zero eigenvalues of the form {±zci, ±ic2, . . . , ±iCm}, Ci, . . . , real, and B has the 
representation 

B = ciBi + C2B2 + . . . + CmBm, (33a) 
BiB^ = B,Bi = 0, 2 ^ J, (33b) 
Bf = -Bi. (33c) 

Lemma [1] is essentially Theorem 2.2 of Gallier and Xu (Q), who also noted the immediate corollary 



= J + XI r[(ci^i + C2B2 + ... + CmB^y 



k=l 



00 _ 
k=l 

^c,B, ^ ^C2B, + , , , + gC„B„ - (m - 1)/ 



= / + X ( sin cjBj + (1 - COS 9)5|) . (34) 
i=i 

Thus, the exponential of any skew symmetric matrix has the form of a sum of Euler-Rodrigues 
terms. 

Lemma 2 The elements of the decomposition of Lemma 1 can he expressed in terms of the matrix 
B and its eigenvalues by 
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Equation (|35|1 may be obtained by starting with 



B{B' + cl) = J2i4-c',)c,B,, (36) 



and iterating until a single Bj remains. Equation fl35l) in combination with flM|) and i?^ = c - ^-B-Bj 
implies that e can be expressed as a polynomial of degree 2m in B, 



i=i k=i S' 



We are now ready to consider exp( 
Lemma 3 The non-zero eigenvalues of the skew symmetric tensor P defined in eq. ^ are 

{i, —i, 2i, —2i, . . . ,ni, —ni}. 



This follows from Zheng and Spence (1231 ) or Lu and Papadopoulos (1l4l). and is discussed in detail 
in Section m P therefore satisfies the characteristic equation ( I2T1) . and has the properties 

P = Pi + 2P2 + . . . + r2P„, (38a) 
PiP,- = P,Pi = 0, z^j, (38b) 
P^ = - Pi . (38c) 

Equation fl38al) is the canonical decomposition of P into orthogonal components each of which 
has properties like the fundamental P in ([3]), although the subspace associated with Pj can be 
multidimensional. The polynomials Pn,k{x) of eq. fll7p follow from fl5B]) . or alternatively, 

We will examine the particular form of the Pj for elasticity tensors of order n = 2 and n = 4 in 
Section [5l For now we note the consequence of Lemmas [2] and [3l 

e^^= I + V [k sin fcg P + (1 - cos k9) P^] fc-^ TT ( 2 {2 ) ■ (^0) 
k=i 3=1 

Together with (13T!) this implies the first expression (I18ap in Theorem [H The alternative identities 
f ll8bp and fll8cp are derived in the Appendix. 

Lemma 4 Q may be expressed in Cayley form 

Q = (1+ S)(I- S)-^ = (I- §)-^(I+ §), (41) 

where the 2n-th order skew symmetric tensor S is 

n 

§ = tanh(| P) = J] tan(fcf ) P^ . (42) 

k=l 
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This follows by inverting (HTj) and using the expression (fT3l) for 



n 

1-2(1+ Q)~^= I- [l + Y^ -(sinA;^ Pfc + (1 -cosA;^)P^)]"\ (43) 



2 

k=l 



and then applying the following identity, 



[ I + E (a. P. + (1 - 6,) pa] " = I + E ( - ^ + (1 - ^) • (44) 



fc=l fe=l 



Note that Lemma H] combined with the obvious result ( H2|) i implies that, for arbitrary 0, 

n 

tanh(</)P) = ^ tan(fc</)) P^. (45) 



fc=i 



This is a far more general statement than the simple partition of (fT2l) . and in fact may be shown 
to be equivalent to eq. (l38ll . 



3.2 Equivalence with the Lu and Papadopoulos formula 



We now show that (140|) agrees with the polynomial of Lu and Papadopoulos (jl4l ). They derived eq. 
(ITTl) directly using Sylvester's interpolation formula (p. 437 of Horn and Johnson (jl3l )). Thus, the 
function /(A) defined by a power series for a matrix A that is diagonalizable can be expressed in 
terms of its distinct eigenvalues Ai, A2, . . . as /(A) = p(A), where p{x) is 

p{x) = j2 /(AO Hx), Hx) = n (^^) . (46) 

1=1 jr = l * ^ 

Li{x) are the Legendre interpolation polynomials and p{x) is the unique polynomial of degree r — 1 
with the property p{\i) = /(Aj). Equation f|TT]) follows from the explicit form of the spectrum of 
P. The right member of ffTTj) can be rewritten in the following form by combining the terms e^*'^^, 

(^0^ ^^J^ [n - k)\ [n + k)\ ^ 

The first term on the right hand side can be written in a form which agrees with fllSap by considering 
the partial fraction expansion of 1/A where 

n 

A{x) = Y[{x^ + . (48) 

k=l 

Thus, as may be checked by comparing residues, for example. 



-=y 



A {n- k)\ {n + k)\ x^ + P 



V 2(-l)^+^ ^ 2(-l)^+^ 

^^{n-k)\{n^k)\ ^ (n-fc)!(r2 + A;)!x2 + P ■ ^ ' 
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The first term in the right member is 1/A(0) = l/(?7,!)^. Multiplying by A(a;)/A(0) implies the 
identity 



(n!)2 Al^ ' ^ ' (n-k)\in + k)\ 

which combined with (1471) allows us recover the form (I18al) . This transformation from the interpo- 
lating polynomial (1111) to the alternative forms in eq. (1181) is rigorous but it does not capture the 
physical basis of the latter. We will find the identity (|50|) useful in Section [61 

4 Eigenvalues of P and application to second order tensors 

The key quantity in the polynomial representation of Q is the set of distinct non-zero eigenvalues 
of P, Lemma [3l We prove this by construction, starting with the case of n = 1. 

4.1 Rotation of vectors, n = 1 

Let {p, q, r} form an orthonormal triad of vectors, then the eigenvalues and eigenvectors of P are 
as follows 

eigenvalue eigenvector 

P (51) 

±i v± = ^(iq±r) 

These may be checked using the properties of the third order alternating tensor. Thus, PijQj = ri, 
PijTj = —Qi, implying Pv-t = ±i\'±. Hence, the spectral representation of P is 

P = iv+v; - iv.vl , (52) 

where * denotes complex conjugate. It also denotes transpose if we view eq. (1521) in vector/matrix 
format. We can also think of this as a tensorial representation in which terms such as v_vl stand 
for the Hermitian dyadic v_ (8> vl, however, for simplicity of notation we do not use dyadic notation 
further but take the view that dyadics are obvious from the context. 

The bivectors (0) v±, together with the axis p, will serve as the building blocks for spectral 
representations of P for n > 1. 

Referring to eq. fl38|) . we see that P = Pi in this case. Equation fl52|) allows us to evaluate 
powers and other functions of P. Thus, in turn, 

p2 = -v+v; - v_vl , (53a) 
p2 + I = pp, (53b) 
P(P' + I) = 0, (53c) 

from which the characteristic equation ([3]) follows. 

4.2 Rotation of second order tensors, n = 2 

In this case the tensors P and Q are fourth order with, see eq. ([9]), 

Pijkl = Pik^jl + SikPjl ■ (54) 
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Equation ([T7|l implies 

Pi = P(p2 + 4)/3, P2 = -P(p2 + l)/6, (55) 
which clearly satisfy the decomposition (I38al) 

P = Pi + 2P2. (56) 

The skew symmetric tensor S of (fT6|) may be expressed 

§ = (I tan I - I tan f ) P + (| tan f - | tan f ) Pl (57) 

Consider the product of the second order tensor (dyad) pv-t with P. Thus, Ppv-t = ±ipv± 
which together with Pv-tp = izv-tp yields 4 eigenvectors, two pairs with eigenvalue i and —i. 
The dyadics v-tV-t are also eigenvectors and have eigenvalues ±2i. The remaining 3 eigenvectors 
of P are null vectors which can be identified as pp, v_|_v_ and v_v+. This completes the 3" = 9 
eigenvalues and eigenvectors, and shows that the non-zero eigenvalues are {z(2), — z(2), 2z, — 2z} 
where the number in parenthesis indicates the multiplicity. The tensors Pi and P2 can be expressed 
explicitly in terms of the eigenvectors, 

Pi = zpv+pv^ + iv+pv^p — ipv_pvl — iv_pvlp, P2 = iv+v+v^v^ — iv_v_vlvl. (58) 

It is straightforward to show that these satisfy the conditions (138|) . 

We note in passing that zero eigenvalues correspond to tensors of order n with rotational or 
transversely isotropic symmetry about the p axis. This demonstrates that there are three second 
order basis tensors for transversely isotropic symmetry. We will discuss this aspect in greater detail 
in Section [6] in the context of elasticity tensors. 

Both P and Q correspond to 3" x 3" or 9 x 9 matrices, which we denote P and Q, with the 
precise form dependent on how we choose to represent second order tensors as 9-vectors. To be 
specific, let Tij be the components of a second order tensor (not symmetric in general), and define 
the 9-vector 

T = (Til, T22, ^33, T23, T31, T12, T32, Ti3, T2i)\ (59) 
The form of P follows from the expansion ([7]) for small 6, or from (Q, as 
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(60) 



This can be written in block matrix form using the X, Y, Z matrices defined in eq. (125|) . 

/O Y Y\ 

P = R-R\ where R = \z X . (61) 

\Z X 0/ 

Note that this partition of the skew symmetric matrix is not unique. 
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4.3 Tensors of arbitrary order n > 2 

The example of second order tensors shows that the three eigenvectors {p, v.,., v_} of the fun- 
damental matrix P define all 3" = 9 eigenvectors of P. This generalizes to arbitrary n > 2 by 
generation of all possible n-tensors formed from the tensor outer product of {p, v+, v_}. For any 
n, consider the product P with the n-tensor 

(62) 

Each of the n terms of P in (Q contributes +iv^ with the result Pw_|_ = niv^. Similarly, defining 
implies Pw_ = —ni w_ , from which it is clear that there are no other eigenvectors with eigenvalues 
±m. Thus, the final component in the decomposition fl38al) is 

P„ = i v+v'^ — i v^v*_ . (63) 

Next consider the n n-th order tensors that are combinations of (n — 1) times v_|_ with a single 
P, e.g. 

u = p v_|_v_|_ . . . v+ . (64) 

n-l 

The action of P is Pv = {n — l)iv, and so there are 2n eigenvectors with eigenvalues ±(n — 1)^ which 
can be combined together to form the component Pn-i, as in eq. (l58|) i for n = 2. Eigenvalues 
(n — 2)z are obtained by combining p twice with v+ (n — 2) times, but also come from n-tensors 
formed by a single v_ with (n — 1) times v^.. Enumeration yields a total of n{n + 1) eigenvectors 
associated with P„_2- By recursion, it is clear that the number of eigenvectors with eigenvalues ki 
equals the number of ways n elements chosen from {—1, 0, 1} sum to —n < k < n, or equivalently, 
the coefficient of x'' in the expansion of (l + a; -|- x~^) , 



(i+^+^-T-i: (J 



. (65) 



The numbers (^)_^ are trinomial coefficients (not to be confused with g-multinomial coefficients ([l| 
with q = 3) and can be expressed (0) 



j=0 



where (") = n!/(j!(n — j)!) are binomial coefficients. The main point for the present purpose is 
that 

k. 



k= 



3", (67) 



2 



which is a consequence of (1651) with x = 1. It is clear from this process that the entire set of 3"" 
eigenvectors has been obtained and there are no other candidates for eigenvalue of P. The precise 
form of the eigenvectors is unimportant because the Euler-Rodrigues type formulas (fT8|) depend 
only on the set of distinct non-zero eigenvalues. However, the eigenvectors do play a role in setting 
up matrix representations of P for specific values of n, and in particular for tensors with symmetries, 
as discussed in Section [51 
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5 Applications to elasticity 

We first demonstrate that the underlying symmetry of second and fourth order tensors in elasticity 
implies symmetries in the rotation tensors which reduces the number of independent elements in 
the latter. 



5.1 Rotation of symmetric second and fourth order tensors 

We consider second and fourth order tensors such as the elastic strain e {n = 2) and the elastic 
stiffness tensor C (n = 4) with elements Eij and Cijki- They transform as eij — > e[j and Cijki — C*^fc;, 

^ij Qijpq^pq ^ ijkl Qijklpqrs Cpqrs ; (68) 

where the rotation tensors follow from the general definition ([8]) as 

Qijpq QipQjq: Qijklpqrs Q ijpqQ klrs • (69) 

The strain and stiffness tensors possess physical symmetries which reduce the number of independent 
elements that need to be considered, 

^ij ~ ^jiy Cijki = Cjikl = Cijlk, Cijki = Cklij- C^o) 

In short, the number of independent elements of e is reduced from 9 to 6, and of C from 81 to 21. 
Accordingly, we may define variants of Qijki and Qijklpqrs with a reduced number of elements that 
reflect the underlying symmetries of elasticity. Thus, the transformation rules can be expressed in 
the alternative forms 

^ij Qijpq^pqi Cijf^i Qijklpqrs Cpqrs i {^^) 

where the elements of the symmetrized 4th and 8th order Q tensors are 

Qijpq ~ 2 {^^ip<i ~^ Qijgp) ' (72a) 

Qijklpqrs ~ 2 {QijpqQ klrs QijrsQklpq) ■ ("^^b) 

Alternatively, they can be expressed in terms of the fundamental Q G 5*0(3), 
_ 1 

Qijpq '^(^QipQjq ~^ QiqQjp)) (73a) 

Q ijklpqrs g (^Q ipQ j qQ krQ Is ~l~ QipQ jqQksQlr ~l~ Q iqQ jpQ krQ Is ~l~ QiqQ jpQ ksQlr 

~\~ QirQjsQkpQlq ~l~ QirQjsQkqQlp ~\~ QisQjrQkpQlq ~\~ QisQ jrQkqQlp) ■ (73b) 

Using the identity (l32l) . we have 

— — 1 - 1 , 

QijpqQ klpq = 2{^ijkl + lijlk) = Ujkl = -^{^ik^jl + ^il^jk), (74) 

which is the fourth order isotropic identity tensor, i.e. the fourth order tensor with the property 
s = Is for all symmetric second order tensors s. Similarly, 

Q ijklabcdQ pqrsabcd ~ 2 {-^ ijPQ-^ klrs + lijrs^klpq) = I ijklpqrs ■ C''^) 
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This is the 8th order isotropic tensor with the property C = IC for all 4th order elasticity tensors. 
In summary, the symmetrized rotation tensors for strain and elastic moduli are orthogonal in the 
sense that Q Q = I, and are therefore elements of 5*0(6) and 50(21), respectively. 
The (symmetrized) rotation tensor for strain displays the following symmetries 

Qijpq Qjipq Qijqp- i'^^) 

Introducing the Voigt indices, which are capital suffices taking the values 1, 2, . . . , 6 according to 

7 = 1,2,3,4,5,6 ^ = 11,22,33,23,31,12, (77) 

then (!76l) implies that the elements Qijpg can be represented by Qjj. Similarly, the elements of the 
8th order tensor Qij^ipq^s can be represented as Qukl^ "which satisfy the symmetries 



Q 



IJKL — "^JIKL 



Q 



IJLK- 



(7J 



The pairs of indices IJ and KL each represent 21 independent values, suggesting the introduction 
of a new type of suffix which ranges from 1, 2, . . . , 21. Thus, 

J = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21 
^ IJ= 11, 22, 33, 23, 31, 12, 44, 55, 66, 14, 25, 36, 34, 15, 26, 24, 35, 16, 56, 64, 45. (79) 

Hence the elements Qijkipqrs can be represented uniquely by the (21)^ elements Qjj- The elements 
of the symmetrized rotations Qjj and Qjj have reduced dimensions, 6 and 21 respectively, but do 
not yet represent the matrix elements of 6- and 21-dimensional tensors. That step is completed 
next, after which we can define the associated skew symmetric generating matrices and return to 
the generalized Euler-Rodrigues formulas for 6- and 21-dimensional matrices. 



5.2 Six dimensional representation 

We now use the isomorphism between second order tensors in three dimensions, such as the strain 
£, and six dimensional vectors according to e — ^ £ with elements ij, I = 1,2, ...,6. Similarly, 
fourth order elasticity tensors in three dimensions are isomorphic with second order positive definite 
symmetric tensors in six dimensions C with elements c/j (jl5l ). Let {e} be the 6- vector with elements 
£/, / = 1, 2, . . . , 6, and [C] the 6x6 Voigt matrix of elastic moduli, i.e. with elements c/j. The 
associated 6-dimensional vector and tensor are 



T{e}, C = T[0]T, where T = diag(l, 1, 1, V2, ^2, V2). 



Explicitly, 















£3 




^33 












V2£3i 









c 



/Cll C12 Ci3 22Ci4 22Ci5 22Ci6\ 

C22 C23 22C24 22C25 22C26 

C33 22C34 22C35 22C36 

2C44 2C45 2C46 

S Y M 2C55 2C56 

V 2C66 / 



(80) 



The \/2 terms ensure that products and norms are preserved, e.g. CijuCijki = tr C*C. 
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The six-dimensional version of the fourth order tensor Q^jf^i is Q G 50(6), introduced by 
Mehrabadi et al. (16), see also (19). It may be defined in the same manner as Q = S[Q]S 
where [Q] is the matrix of Voigt elements. Thus, QQ* = Q*Q = I, where I = diag(l, 1, 1, 1, 1, 1). 
It can be expressed Q(p,^) = exp(6'P), and hence is given by the n = 2 Euler-Rodrigues formula 



22l) for the skew symmetric P(p) G so(6). 



/ 











V2P2 


-V2P3\ 











-V2pi 





V2p3 











V2pi 


-V2p2 








V2pi 


-V2pi 







-P2 


-V2p2 





V2p2 


-P3 





Pl 


\ V2P3 


-V2p3 





P2 


-Pl 


/ 



(82) 



or in terms of the block matrices defined in 



R-R* 



R 



V2Y 
V2Z X 



(83) 



It is useful to compare P with the 9x9 matrix for rotation of general second order tensors, eqs. (l60l) 
and (EH) . The reduced dimensions and the a/2 terms are a consequence of the underlying symmetry 
of the tensors that are being rotated. Vectors and tensors transform as e e and C ^ C where 



e' = Qe, C' = QCQ*. 
The eigenvalues and orthonormal eigenvectors of P are as follows 

dyadic 



(84) 



eigenvalue 




±i 

±2i 



eigenvector 
i 

Pd 
u± 
v± 



|(PP 



ID 

^(pv± + v±p) 
v±v± 



(85) 



where v-t are defined in (I5T!) and 

/pI- 



i/Vs 



V / 



Pd 



pI 



i\ 

1 

3 
1 



V2p2P3 

V^psPi 

W2P1P2J 



/ y/2piv±^i \ 

V^P2V±,2 
V2p3V±,3 
P2V±,3 + P3V±,2 
P3V±,1 +PlV±,3 
\PlV±,2+P2V±,lJ 



v± = 



"±,1 

^2 

As 



\/^V±,2V±,3 
VV2f±,lt'±,2/ 



• (86) 



The right column in ( l85l) shows the second order tensors corresponding to the eigenvectors in dyadic 
form. The eigenvectors are orthonormal, i.e. of unit magnitude and mutually orthogonal. The unit 
6— vectors i and correspond to the hydrostatic and deviatoric parts of pp, respectively. The 
double multiplicity of the zero eigenvalue means that any 6— vector of the form ap^ + 6i is a null 

vector of P, including y fpd + which is the six- vector for the dyad pp. 
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The canonical decomposition of the generating matrix P e so{6) is therefore 

P = Pi + 2P2, Pi = m+u;; -m_ul, P2 = iv+v; - iv_vl . (87) 
The associated fourth order tensor P, which is the symmetrized version of P in eq. (1561) . is 



P=Pi + 2P2, Pi = J] ^(pv± + v±p)(pv^ + vip), P2 = J] ±zv±v±v;vi. 



The representation fISTI) allows us to compute powers of P, and using the spectral decomposition 
of the identity 

I = ii* + PdPd + u+u* + u_ul + v+v;^ + v_vl , (89) 



the characteristic equation (jl6l ) follows: 

P(PVI)(P^ + 4I) = 0. (90) 

5.3 21 dimensional representation 

Vectors and matrices in 21-dimensions are denoted by a tilde. Thus, the 21-dimensional vector of 
elastic moduli is c with elements cj, X = 1,2,..., 21. Similarly, vectors corresponding to the fourth 
order tensors ££ and ^(A ® B + B (g) A) are £ and v, that is 

^th Qj.(^Qj.^ n=3 2"''^ order, n=6 vector, n=21 

i(A®B + B®A) |(ab* + l3a*) v 
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The 21 elements follow from the indexing scheme (179|) and are defined by 







( 1 




( 1 




( \ 




( \ 


C2 




C22 




C22 




^22 




el 


C3 




C33 




C33 




^33 




el 


C4 




\/2C23 




\/2c23 




■v/2e22£33 




\/2e2e3 


C5 




\^Ci3 




\^Ci3 




V2e 11^33 






C6 




\/2ci2 




\^Ci2 




\/2eiie22 




\/2eie2 


C7 




2C44 




C44 




^^23 




^1 


C8 




2C55 




C55 




2p2 
^^13 




^1 


eg 




2C66 




C66 




^^12 




il 


ClO 




2ci4 




\/2ci4 




2eiie23 






Cll 


= 


2C25 


= 


V2C25 


£ = 


2e22ei3 


= 




Cl2 




2C36 








2e33£l2 




\/2e3e6 


CIS 




2C34 




\/2C34 




2^33 £23 




\/2e3e4 


Ci4 




2ci5 




\/2ci5 




2eiiei3 




\/2eie5 


Cl5 




2C26 




\^C26 




2e22ei2 




\/2e2e6 


Cl6 




2C24 




\/2c24 




2e22e23 




\/2e2e4 


Cl7 




2C35 




\/2C35 




2e33ei3 




\/2e3e5 


Cl8 




2C16 




\/2ci6 




2eiiei2 




\/2eie6 


Cl9 




2\/2c56 




\/2c56 




2\/2ei3ei2 




\/2e5e6 


C20 




2\/2c46 




\/2c46 




2i/2e23£i2 




\/2e6e4 






l^2V2C45yi 




i^\/2c45 / 




i^2^/2e23ei3/ 




V2e4e5 / 



Alternatively, 



02^2 

03^3 
:^(a2^3 + 03^2) 
^(0361 + 0163) 

^(0162 + a2&i) 
0464 
0565 



-^(0164 + 0461 
^(0265 + 05^2 
^(0366 + a6&3 
^(0364 + 0463 
^(0165 + 0561 
^(0266 + 0662 
^(0264 + 0462 
-^(0365 + 0563 
^(ai^e + aefti 
^(0566 + 0665 
-^(0664 + 04^6 

+ 05^4) J 
(92) 



(93) 



c = T{c}, etc., 

where {c} is the 21x1 array with elements cj, and T is the the 21 x 21 diagonal with block structure 

T = diag (I V2I 21 21 21 21 2^21) . (94) 
For instance, the elastic energy density W = \Cijkieijeki is given by the vector inner product 

1 



W 



where e is the 21- vector associated with strain 



-t~ 

£ C 



(95) 



5.3.1 The 21-dimensional rotation matrix 

Vectors in 21-dimensions transform as c — c' = Qc where Q gS0(21) satisfies QQ* = Q*Q = I, 
and I is the identity, Ijj = Sjj. The rotation matrix Q can be expressed 

„9P 



Q(p, 



(96) 



where P(p) G so(21) follows from the general definition (Q applied to elasticity tensors, i.e. (173bp . 

Thus, _ _ _ _ _ _ 

P = T[P]T, Q = T[g]T (97) 
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where [P], [Q] are the 21 x 21 array with elements Pjj, Qjj- We take a slightly different approach, 
and derive the elements of P using the expansion of Q for small 9 directly. Define the rotational 
derivative of the moduli as 



Ajkl 



dc'ijki 



d9 



Using fl8^ 9 for instance, gives 



C = PC + CP* 



from which we obtain 

Cll = 

C22 = 
C33 = 
C23 = 

Cl3 = 
C12 = 
C44 = 

C55 = 

Cm = 

Ci4 = 
C25 = 
C36 = 
C34 = 
Cl5 = 
C26 = 
C24 = 
C35 = 
C16 = 
C56 = 
C46 = 
C45 = 



4ci5P2 - 4ci6P3 , 
-4C24Pl + 4C26P3 , 
4C34P1 - 4C35P2 , 

2(C24 - C34)pi - 2025^*2 + 2C36P3 , 

2CuPi + 2(C35 - Ci5)j92 " ^C^ePs , 

-2Ci4Pl + 2C25P2 + 2(ci6 - C26)P3 , 

2(C24 - Cu)Pl - 2C46P2 + '^C^bPs , 

2C56P1 + 2(C35 - Ci5)p2 " ^C^bPs , 

-2C56P1 + 2C46P2 + 2(ci6 " C26)P3 , 

(Cl2 - Ci3)j9i - (Ci6 - 2C45)P2 + (Cl5 " 2C46)p3 , 

(C26 - 2c45)pi + (C23 - 012)^2 " (C24 - 2C56)p3 , 

-(C35 - 2C46)P1 + (C34 - 2C56)P2 + (Cl3 - £23)^3 , 

-(C33 - C23 - 2c44)pi - (C36 + 2C45)P2 + C35P3 , 

CwPl - (Cii - Ci3 - 2C55)P2 - (Cl4 + 2c56)p3 , 

-(C25 + 2C46)P1 + C24P2 - (C22 " C12 - 2cqq)P3 , 

(C22 - C23 - 2C44)pi - C26P2 + (C25 + 2C46)p3 , 

(C36 + 2C45)P1 + (C33 - Ci3 - 2C55)P2 " C34P3 , 

-ClsPl + (Ci4 + 2C56)P2 + (Cll - C12 - 2C66)P3 , 

(C66 - C55)Pl + (C36 " C16 + 045)^2 " (C25 " C15 + 046)^3 , 

-(C36 - C26 + C45)pi + (C44 - C6e)p2 + (cu - C24 + £56)^3 

(C25 - C35 + C46)pi - (Ci4 - C34 + C56)p2 + (C55 " £44)^3 • 



(98) 



(99) 



(100a) 
(100b) 
(100c) 
(lOOd) 
(lOOe) 
(lOOf) 
(lOOg) 
(lOOh) 
(lOOi) 
(lOOj) 
(100k) 
(1001) 
(100m) 
(lOOn) 
(lOOo) 
(lOOp) 
(lOOq) 
(lOOr) 
(100s) 

(loot) 

(lOOu) 



The elements of P can be read off from this and are given in Table 1. The 21 x 21 format can be 
simplified to a 7 x 7 array of block elements, comprising combinations of the 3x3 zero matrix 



and the 3x3 matrices defined in 


m 


Thus, 














(0 











2Y 





\ 













-V2Y 





v^N 






















2N 


-V2Y 


P = R - R*, where R = 





-V2Z 








X 





-V2X 







v^N 


2N 








X 







2Z 








X 

















-V2Z 


-V2X 


V2X 








(101) 
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Equation flM|) shows that the rotation matrices in three, six and 21 dimensions can be simphfied 
using the fundamental 3x3 matrices X, Y, Z and N. The partition of a skew symmetric matrix 
in this way is not unique, but it simphfies the calculation and numerical implementation. These 
matrices satisfy several identities, 

XX* + YY* + ZZ* = I , (102a) 

X*X = YY*, Y*Y = ZZ*, Z*Z = XX*, (102b) 

X*Y = YZ*, Y*Z = ZX*, Z*X = XY*, (102c) 

||N|| = ||X|| = ||Y|| = ||Z|| = 1, (102d) 



where the norm is = yirvFu . 

We find that the eigenvalues of P are 0(5), ^(3), — i(3), 2i(3), — 2z(3), 3i, — 3i, 4i, —Ai, where 
the number in parenthesis is the multiplicity. The associated orthonormal eigenvectors of P are 

eigenvalue eigenvector 6-dyadic 
voa ii* 

Vo6 fPdPd 

yo, ^(ip*i + pj*) 

Vod ^(u_u+ + u+u*„) 

Voe ^(V_V^^. + V+V*L) 

±« via± ^(iu*i + u±i*) (103) 

±^ yi6± ^(PdU*± + u±p^) 

±^ Vie± -lj(v±U*^ + U^Vy 

±2i V26± ^(iv^ + v±i*) 

±2^ y2,± ^(PdV±+v±p*,) 

±32 V3± ;7|(u±v^ + v±u^) 

Apart from the last two, which have unit multiplicity, these are not unique. For instance, different 
combinations of the eigenvectors with multiplicity greater than one could be used instead. We also 
note from their definitions in eqs. fl5T]) and fl85|) that 

v; = -v^, = -u^, = v^. (104) 

Hence, the five null vectors of P in fll03p are real. P can be expressed in terms of the remaining 16 
eigenvectors as 

P = Pi + 2P2 + 3P3 + 4P4, (105) 

where 

P, = Y.^'ll ^j^±^U^ J = 1, 2; P,- = v,±v*i, J = 3, 4. (106) 

it x=a,b,c 

The identity is 

1= Yl ^Ox^Ox + ^iYl (Vlx±vt^± + V2^±V;^±) + V3±V*± + V4±V*±) . (107) 

x=a,b,c,d,e 
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It is straightforward to demonstrate using these expressions that P satisfies 

P (P2 + 1) (P2 + 4) (P2 + 9) (P2 + 16) = , (108) 
which is the characteristic equation for 2?7,-th order P tensors, eq. (12T|) . 

6 Projection onto TI symmetry and Cartan decomposition 

It has been mentioned in passing that null vectors of P and its particular realizations for elasticity 
tensors correspond to tensors with transversely isotropic symmetry. These are defined as tensors 
invariant under the action of the group 5*0 (2) associated with rotation about the axis p. We 
conclude by making this connection more specific, and relating the theory for the rotation of tensors 
to the Cartan decomposition, which is introduced below. 

6.1 Cartan decomposition 

Referring to Theorem[31 the expansion of the projector defined in eq. f l27l) as an even polynomial 
of degree 2n in P, i.e. of degree n in P^, follows from eqs. fl39p and fl50p and the identity P^, = 
^-ipp^_ Note that the Mfc tensors partition the identity 

I = Mo + Ml + . . . + M„. (109) 

Starting with Mq we have 

PMo = MqP = 0, (110) 

indicating that Mq is the linear projection operator onto the null space of P. In order to examine 
the remaining n projectors, note that the rotation can be expressed 

n 

Q= Mo + ^ (cosj^Mj + sinj^ P^) . (Ill) 
i=i 

The action of Q on an n-th order tensor T, the rotation of T, can therefore be described as 

n 

QT= To + ^ {cosjOTj + smjeRj) , (112) 

where Tq and Tj, Mj, j = 1, 2, . . . , n are defined by (!26l) . We note that may be expressed 

Rj = Pj = - P Tj, J = 1, 2, ... , n. (113) 

Also, To is the component of T in the null space of P, i.e. the transversely isotropic part of T. 
Equation ( 1112^ show that Tq and the pairs Tj, Mj, j = l,2,. . . , n, rotate separately, forming distinct 
subspaces, and hence proving Theorem [3l 

6.2 Application to elasticity tensors 

The projection operators are now illustrated with particular application to tensors of low order and 
elasticity tensors, starting with the simplest case: rotation of vectors. 
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6.2.1 Vectors, n=l 

If n = 1 , we have T = P of eq. ([T]) , and hence 

Mo = pp, Ml = I - pp. 
Also, T ^ t, an arbitrary vector, and the decomposition flll2p and fl26l) is 

Qt = (t ■ p)p + cos 6* [t - (t ■ p)p] + sin 61 p X t. 



(114) 



;ii5) 



Thus, we identify to = Mot as the component of the vector parallel to the axis of rotation, ti = Mit 
as the orthogonal complement, and ri = Pit is the rotation of ti about the axis by 7r/2. 



6.2.2 Second order tensors, n=2 

For n = 2, we have Pi, P2 given in (!55|) and 

Mo = (p2 + l)(p2 + 4)/4, Mi = -P^(p2 
The dimensions of the subspaces are Tq, 3; Ti, 2, ; 



4)/3, 



2; T,, 1, 



p2(p2 + 1)/12. 
1. 



;ii6) 



For symmetric second order tensors the projectors Mq, Mi and M2 have the same form as in 
(11161) in terms of P, and may also be expressed in terms of the fundamental vectors introduced 
earlier. Thus, using (!87|) and (!89|l . 



Mn 



4 



I)(p2 + 4I) = ii* + PdM 



Ml 



Mo 



V V 



(117) 



Using these and Pi, P2 from (!87|) . it follows that the subspace associated with Mq is two dimen- 
sional, while Mi^\ Mp"*, and Mj^^ have dimension one. This Cartan decomposition of the 



6-dimensional space Sym agrees with a different approach by Forte and Vianello ((8|). Alternatively, 
using (11171) and the dyadics in ( l85l) . it follows that 



Mo = pppp + ^(I - pp)(I - pp), 

1 1 
Ml = -(pq + qp)(pq + qp) + -(pr + rp)(pr + rp), 

M2 = ^(qq rr)(qq - rr) + ^(qr + rq)(qr + rq). 

The action of Mq on a second order symmetric tensor is 

MoT = To = 3||pp + rx(I-pp), 



(118a) 
(118b) 
(118c) 

(119) 



where Ty = T : pp and T± = i(trT — T : pp). It is clear that Tq is unchanged under rotation 
about p. Similarly, MiT transforms as a vector under rotation, hence the suffix 1, while M2T 
transforms as a second order tensor with elements confined to the plane orthogonal to the axis of 
rotation. 

Let p = e3, then using (I117l) i or the simpler (I118ap yields 



Mo(e3)T 



/I 0\ 

i 

1 





\0 0/ 



/(fi+f2)/2\ 

(ri + f2)/2 

T3 




V / 



(120) 
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The projection To can be identified as the part of T which displays hexagonal (or transversely 
isotropic) symmetry about the axis 63. To be more precise, if we denote the projection MqT as 
Tnex; then Thox is invariant under the action of the symmetry group of rotations about p. It 
may also be shown that of all tensors with hexagonal symmetry Thcx is the closest to T using a 
Euclidean norm for distance. Similarly, 

ti = Mit = (0, 0, 0, f4, % 0)*, % = Mst = (l(fi - fs), ^(fs - fi), 0, 0, 0, %)\ (121) 

and using eqs. (l82l) and fll2ip . 

Ri = (0, 0, 0, fs, -Ti, 0)*, R2 = ( - -^Te, -^fg, 0, 0, 0, -^(fi - f^))'. (122) 

It can be checked that under rotation, we have 

QTo = to, QTj- = cos j^tj + sin j^Rj, QR^ = cos j^R^ - sin j^f j-, j = 1, 2. (123) 

Thus, the singleton to and pairs (ti, Ri) and (t2, R2) form closed groups under the action of the 
rotation. 



6.2.3 Elastic moduli, n=4 



The 5-dimensional null space of P is equal to the set of base elasticity tensors for transverse isotropy 
The projector follows from eqs. fll06p and fll07p as 



Mo 



(4!)^ 



;{P^ + 1) (p2 + 4) (p2 + 9) (p2 + 16) 



(124) 



x=a,b,c,d,e 



This is equivalent to the 8-th order projection tensor Phbx of eq. (68c) of Moakher and Norris (1171). 
who use the base tensors suggested by Walpole (1221 ). The latter representation is useful because 
Walpole's tensors have a nice algebraic structure, making tensor products simple. Alternatively, in 
the spirit of the 21 x 21 representation of elastic moduli, (11241) 1 and (11011) yield 
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(125) 
w for 



The projection Mq of (1125p is identical to the matrix derived by Browaeys and Chevrot 
projection of elastic moduli onto hexagonal symmetry. Specifically^the 9x9 matrix Mhcx is exactly 
M*-^-* of Browaeys and Chevrot (see Appendix A of (p)). Hence, Mq of (1125p provides an algorithm 
for projection onto hexagonal symmetry for arbitrary axis p. 

The subspace of To is 5-dimensional, while the subspaces (Ti, Ri), (T2, R2), (T3, R3) and (T4, i?4) 
have dimensions 6, 6, 2 and 2, respectively. These may be evaluated for arbitrary axis of rotation 
p using the prescription in eqs. ( 128cl) and (11131) . 
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A. Appendix: A recursion 



Equation fl371) provides an expansion in terms of basic trigonometric functions. An alternative 
recursive procedure is derived here. Consider the expression fl37|) written 



e'' = fmiB), (A.l) 

where the function fm{B) is a finite series defined by the m distinct pairs of non-zero eigenvalues 
of B. Suppose the set of eigenvalues is augmented by one more pair, ±icm+i, with the others 
unchanged. Then the sums and products in (1371) change to reflect the new eigenvalues, 

m+l m+1 „2 I 2 

= fm+i{B) =I + J2 si^c^- B + cose,) B^] cf J] { 2 _ S ) ■ (^-2) 

j=i k=i '^1' 

Mi 

The effect of the two additional eigenvalues can be isolated using the identity 
to give 

This defines the sequence, starting with /q = /. 

Applying this to the rotation exp(^P), and noting that the eigenvalues of the skew symmetric 
tensor P correspond to cj = jO, gives 

JSmj^X + (1 - COSJ^jX \ -p-r 2 , ,2n / a c;^ 

2^ ~i (A.5) 

1=1 

The sum over trigonometric functions may be simplified using formulas for sinn^ and cosnO in 
terms of cos6/2 powers of sin^/2 in Section 1.33 of Gradshteyn et al. (jiol ). Thus, 

^ ^^+1 - (2n+l)! '°'2' ^"^-^^ 



(=1 



n - j') 

and hence, 

\2n+l 



y 1-cosjg _ (2sinf)^"+^ sinf 

^ n+i ~ (2n + l)! n + 1 ' ^ ' ^ 



P2n+2{e, X) = P2n{e, x) + \^ ( COS - X + (n + 1)-' siu -x') ^[(x^ + k') , (A.8) 

(2n + l)! '2 2 ^^Jj 

from which the expression ( 118bl) follows. 
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Table 1. The matrix P(p) G so(21) for 4th order elasticity tensors represented as 21-vectors. The rotation Q = exp(^P) is 

given by the 8-th order polynomial Q = PgiO, P) of eq. fpUj) . 



